Tunable fishnet metamaterials infiltrated by liquid crystals 
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We analyze numerically the optical response and effective macroscopic parameters of fishnet metamaterials 
infiltrated with a nematic liquid crystal. We show that even a small amount of liquid crystal can provide 
tuning of the structures due to reorientation of the liquid crystal director. This enables switchable optical 
metamaterials, where the refractive index can be switched from positive to negative by an external field. This 
tuning is primarily determined by the shift of the cut-off wavelength of the holes, with only a small influence 
due to the change in plasmon dispersion. 



Metamaterials have attracted a great attention due to 
their unusual electromagnetic properties not available in 
natur^. One of the important metamaterial designs that 
is suitable for scaling to optical frequencies is the fishnet 
metal-dielectric-metal structur^^H^. 

Due to the strong local field enhancement within the 
metamaterials, there emerge new opportunities for ex- 
ternal control of their properties or changing their op- 
erating bandwidth. There are several means to achieve 
such tunability, including structural^, temperatur d^ * -^^*, 
electro-opt iJ^^J^, magneto-optiJ^, and nonlinear tun- 
^j^gi4"i5'^ However, so far there are no works on tuning 
the optical properties of realistic structures with negative 
index, such as the fishnet structure. 

In this Letter, we numerically analyze the tunability 
of optical fishnet metamaterials exhibiting negative re- 
fractive index. We study a structure with the holes in- 
filtrated by liquid crystals and demonstrate that even a 
small amount of liquid crystal can allow for substantial 
tunability of such structures, including reversal of the 
sign of the refractive index. This is achieved by employ- 
ing the reorientation of the liquid crystal director using 
an external field, as well as by implementing the temper- 
ature dependence of the liquid crystal refractive index. 

Liquid-crystal tunability of metamaterials has been 
suggested as a path to achieve tunable negative- index 
metamaterials^^. Subsequent studies have explored this 
idea in quasi-two-dimensional nanostrip geometrjE'^ 
and also for a wire medium^^. Recent experiments 
demonstrated thermal tunability of optical metamateri- 
als, however the realized structures only experienced a 
negative magnetic response, and did not show a negative 
refractive index^. Therefore, here we present a study on 
tuning of optical fishnet structures as an important step 
in achieving tunable negative-index optical metamateri- 
als. 

Following our recent experiments^, we consider a 
trilayer metal-dielectric-metal structure patterned with 
a two-dimensional square lattice of rectangular holes 
through all three layers. The structure is deposited on a 
glass substrate, as shown in Fig. [l] In our simulations, 
we use realistic material parameters as gold (Au) and 
zinc-oxide (ZnO) for the metal and the dielectric layers, 
respectively. We fix the total thickness of the structure, 
h = 2hm -\-hd = 120 nm, where hm is the thickness of the 
gold and hd is the thickness of the ZnO layers. We use ge- 




FIG. 1. (color online) Trilayer fishnet metamaterial on a 
glass substrate. Left: Schematic of the layers. Right: Four 
unit cells of the structure, km is the thickness of metal (Au), 
hd is the thickness of dielectric (ZnO). Both substrate and 
cladding are semi-infinite. Light is polarized along x axis. 

ometric parameters of the square unit cell of a = 200 nm, 
b = 350 nm, and c = 420 nm (see Fig. [T]). 

We employ a finite-difference time-domain (FDTD) 
numerical method performed with commercial software 
(RSoft). For simplicity, in our simulations we assume 
that the glass substrate and surrounding air are semi- 
infinite with refractive indices ni = 1.48 and na = 1, re- 
spectively. The refractive index of ZnO is taken nznO — 
1.66^^ and its weak dispersion in the spectral range of 
1 — 2 /im (maximal deviation from 1.66 is 1.1%) allows for 
a good approximation by a constant value. For the dielec- 
tric permittivity of gold we use six-term Drude-Lorentz 
expression from the RSoft material librar}!^. A detailed 
description of the implementation of the model in RSoft 
is given in Ref.^. In our simulations, we use a nonuni- 
form spatial grid with size varying from 2nm near the 
metal-dielectric interfaces to 20 nm in free space. The 
time step (in the units of d) is 8 x 10~^ /im. In order to 
simulate a periodic structure in the (x, y) plane we use 
periodic boundary conditions on the planes orthogonal 
to X and y and perfectly matched layer boundary condi- 
tions on the sides orthogonal to z. The structure is ex- 
cited from the air by a short pulse with a broad spectrum 
and plane wave front. Several time monitors are used in 
the simulation to record the electric field components of 
the reflected and transmitted waves. The monitors are 
located 4.4 /im from the structure, so only propagating 
waves are recorded. We take a Fourier transform of the 
response and normalize it to the incident signal, compen- 
sating for the phase shift through air and the substrate. 
The obtained transmission t and reflection r which are 
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FIG. 2. (color online) (a, c) Real part of the effective refrac- 
tive index of the fishnet structure as a function of ZnO layer 
thickness, (b, d) Transmission as a function of ZnO layer 
thickness. Top and bottom row correspond to a director of 
the liquid crystal aligned along z or x axes, respectively. 

referenced to the surfaces of the metal fishnet. 

From the complex coefficients of transmissions and re- 
flection, we can extract the effective refractive index n 
and the impedance z of the metamaterial structure using 
the Fresnel formulae^. For an equivalent isotropic homo- 
geneous slab of thickness h surrounded by semi-infinite 
media with refractive indices ni and ns and normal inci- 
dence, the expressions are 
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FIG. 3. (color online) (a) The real (solid line) and imaginary 
(dotted line) parts of the effective refractive index, (b) Trans- 
mission (solid line) and reflection (dotted line). The thickness 
of ZnO layer hd — 25 nm. (c) Transmission (solid line) and 
reflection (dotted line) of the respected hole array without 
ZnO layer. Black lines correspond to a liquid crystal director 
oriented along z and red lines to a director along x. 
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where m is an integer number. This result agrees with 
the previously derived formulae in the limit ni = ns = 1, 
and the same considerations apply to the choice of signs 
and the branch index tttP^I. We base our choice of branch 
on the requirement that n converges in the low frequency 
limit, and our curves are found to agree with previous 
results using similar geometrical parameters'^ 

We then apply this extraction approach to determine 
the change in effective index of the metamaterial when 
the orientation of the director d is changed by 90°. We 
are most interested in the case when such change of the 
director will switch the refractive index from positive to 
negative. Therefore, we seek the value of the dielectric 
layer thickness which maximizes this effect. We perform 
a number of simulations varying hd and using two values 
of the refractive index of the liquid crystal: = 1.5 for 
ordinary waves when the director d\\ y or d\\ z (Fig. [l] 
left); and ng = 1.7 for extraordinary waves when d || x. 
These values correspond to E7 liquid crystal at 25° O^. 



We note that in the spectral range of 1 — 2 jiia the maxi- 
mum deviation of the liquid crystal refractive index from 
these nominal values is 0.25% for and 0.35% for ng. 
For simplicity, we also simulate the liquid crystal as a ho- 
mogeneous medium with isotropic refractive index. This 
is a reasonable approximation given that the electric field 
within the holes is mostly oriented in the same direction 
as the incident field^. 

The extracted effective refractive index and energy 
transmission coefficient T = n^/ni\t\^ are presented in 
Fig.|2]as a function of hd (varied in 5 nm steps). The neg- 
ative effective index in the fishnet structures is associated 
with the excitation of plasmons between the two metal 
layers'^, and the absorption due to these modes can be 
seen as transmission dips in Fig. |2|b, d). The mode ex- 
cited at longer wavelengths leads to a negative refractive 
index for a range of thicknesses hd- There is also a weakly 
excited mode at shorter wavelengths which can also yield 
to a negative index if other structural parameters are se- 
lected appropriately. Most importantly, we see that the 
effective refractive index switches from negative to pos- 
itive value for ZnO thickness of around 25 nm (shown 



3 



— t=15°C d||z 


-t=15°C d||x 


— t=55°C d||z 


— t=55°C d||x 




?^(|Lim) }i(|Linn) 

FIG. 4. (color online) (a) The real (solid line) and imaginary 
(dotted line) parts of the effective refractive index, (b) Trans- 
mission (solid line) and reflection (dotted line). The thickness 
of ZnO layer is hd — 25 nm. 

by the dashed line). This effect takes place in a narrow 
thickness range of only ~ 5 nm, however this is feasible to 
fabricate with contemporary sputtering techniques that 
can have a precision of 1 nm. 

To understand the effects leading to switching of the 
effective index for 25 nm ZnO thickness, in Fig.|3]we plot 
the extracted refractive index, energy transmission and 
reflection for both orientations of the liquid crystal direc- 
tor. Due to the strong spatial dispersion of the fishnet 
structure^^, we avoid using the local effective parameters 
e and /i, since they serve only to reconstruct n and z and 
lack physical meaning of their own. Instead, we recall 
that the physical origin of the negative index is a gap 
plasmonic mode, with a strong magnetic response, ly- 
ing in the region of plasma-like electrical response of the 
structurd^. In fact, this effective plasma response ex- 
ists above a cut-off wavelength that is close to the cut-off 
wavelength of the fundamental mode of the hole array. 

In order to isolate the influence of the hole modes 
on the response, in Fig. [sj^c) we show the calculated 
transmission through a structure where the ZnO layer 
is replaced by gold, thus suppressing the gap plasmonic 
modes. It can be seen that the transmission peak and 
reflection dip in Fig. [sj^b) at around 1.3 /im correspond 
to the cut-off resonance of the hole array, which changes 
drastically with the change of liquid crystal director ori- 
entation. On the other hand, the transmission dip in 
Fig. [sj^b) that corresponds to the plasmonic dispersion 
curve visible in Fig.[2j barely shifts, since these modes are 
expected to be largely confined between the metal layers. 
Thus, when tilc increases, the cut-off wavelength shifts 
further to the infra-red by approximately 170 nm and the 
hybridization between the hole and plasmon modes which 
led to the negative index is disturbed. 

We also study the effect of temperature on the infil- 
trated metamaterials due to the strong temperature sen- 
sitivity of liquid crystal refractive index, especially near 
the clearing temperature (58° C for E7). This temper- 
ature dependence could potentially disturb the reorien- 
tation based tuning of the metamaterials. In Fig. [4] we 
display the extracted effective refractive index, energy 
transmission and reflection at 15°C {rio = 1.5, rie = 1.71) 



and at 55°C (n^ = 1.52, Ue = 1.65)P for both orienta- 
tions of the liquid crystal director. It can be clearly seen 
that the negative index region undergoes almost no shift 
with temperature, and that at both temperatures we are 
able to switch the index between positive and negative, 
thus the reorientational tuning is robust to temperature. 

In conclusion, we have demonstrated numerically the 
tunability of optical fishnet metamaterials infiltrated by 
nematic liquid crystals. We have found that by tuning 
of the director orientation of an anisotropic liquid crystal 
we can reverse the sign of the effective refractive index 
of the fishnet metamaterials. We have shown that this 
tuning is dominated by a shift in the cut-off frequency of 
the hole modes, and that the negative index is robust to 
temperature variations. We believe that our results will 
facilitate further experimental studies of tunable optical 
metamaterials. 
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